Microscopic Approach to Shear Viscosities in Superfluid Gases: Prom BCS to BEC 
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We compute the shear viscosity, n, at general temperatures T, in a BCS-BEC crossover scheme 
which is demonstrably consistent with conservation laws. The study of ij is important because it 
constrains microscopic theories by revealing the excitation spectra. The onset of a normal state 
pairing gap and the contribution from pair degrees of freedom imply that r\ at low T becomes small, 
rather than exhibiting the upturn predicted by most others. Using the local density approximation, 
we find quite reasonable agreement with just-published experiments. 
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The ultracold Fermi gases near the unitary limit are 
thought to be related to quark-gluon plasmas Much 
attention has focused on particle-physics-based calcula- 
tions of the anomalously low shear viscosity, rj @, 0|. 
However, considerable insight on the thermodynamics 
[H and various spectroscopic studies @, @ has also been 
obtained via a condensed matter perspective. This pa- 
per belongs to the second school in which BCS theory is 
extended to accomodate arbitrarily strong interactions. 
We apply BCS to Bose Einstein condensation (BEC) 
crossover theory [B| to compute rj demonstrating consis- 
tency with central sum rules Q and conservation laws. 

The shear viscosity is a powerful probe for testing mi- 
croscopic theories, because it reflects the normal fluid 
component. As a result it is extremely sensitive to the 
nature of the excitation spectrum. In a low T normal 
Fermi liquid phase with scattering lifetime 7 -1 and ef- 
fective mass to*, rj = ^nvpj^ 1 !?!* . More generally, one 
can think of rj as characterized by the effective number 
of the normal excitations (n — > n e ff(T)) as well as their 
lifetime which we emphasize here is a many body effect. 
Crucial is an understanding of how n e tf depends on T. 
BCS-based approaches have addressed the behavior of 
the viscosity in a notable class of fermionic superfiuids- 
helium-3 @, Q. Here experiments [l(J indicate that 77 
drops off rapidly to zero in the superfluid phase, reflect- 
ing the suppression of fermionic excitations at low T. 
In the helium-4 counterpart, the single particle bosonic 
excitations couple to the collective (Nambu-Goldstone) 
modes, leading to an upturn [11] in 77 at low T, which has 
also been predicted (but not seen) for the atomic Fermi 
supcrfluids [3]. In BCS-based superfiuids, we stress that 
Nambu-Goldstone boson effects do not naturally enter 
into the transverse transport properties, such as rj. 

The Fermi superfiuids involve [l2| inter-dependent 
fermionic as well as bosonic (or pair) excitations, here de- 
riving from stronger-than-BCS attraction. In past litera- 
ture there has been a focus on either one [HI or the other 



3], but not both. Here we use a Kubo-based formalism 
which readily accomodates the simultaneous bosonic and 
fermionic contributions and thereby addresses n e f f quite 
accurately while the alternative Boltzmann or kinetic 
theory-based approaches do not naturally incorporate 
these multiple statistical effects. The Kubo approach 
includes scattering processes via the lifetimes [14| which 
appear in the various Green's functions, while Boltzmann 
schemes treat lifetimes via collision integrals. However, 
because the physics of this dissipation is principally as- 
sociated with the many body processes of boson-fermion 
inter-conversion, it can be satisfactorily addressed only 
in theories which treat the mixed statistics. In this way 
it appears that Boltzmann based schemes may turn out 
to be inadequate, particularly at low T . 

Our central conclusion here is that both the effects of 
a fermionic gap (with onset temperature T* > T c ) and 
the non-condensed pairs act in concert to reduce n e ff 
and thus lower the shear viscosity at all T < T*. When 
quantitatively compared with very recent shear viscosity 
experiments [l| (we independently infer [l5| an estimated 
lifetime from radio frequency data) the agreement is rea- 
sonable. Using previous thermodynamical calculations 
dj of the trap energy E and entropy density s, a plot of 
the trap-integrated rj decreases with decreasing E, while 
rj/ s appears to be roughly constant at low E, but above 
the universal quantum limit 

We compute viscosities using the more well controlled 
current-current correlation functions 0] , XJJ rather than 
stress-tensor correlation functions, via 

n = —to 2 lim lim — Ix3XYt(oJ,ci) 

w-s-o q->o q z 

which is importantly constrained by the sum rule 
du i ImxT(w,qK n n (T) 
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The transverse susceptibility \t = (ELi Xjj ~ Xl)/2 
with the longitudinal xl = q • XJJ ' Q- Here n n {T) cor- 
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responds to the number of particles in the normal fluid 
and n n (T) — > n above T c , where n is the total particle 
number. Because dissipative transport in BCS-BEC the- 
ory is complex, among the most persuasive checks on a 
proper characterization of the "normal fluid" is consis- 
tency with sum rules [25j. Equation ([T]) is more general 
and fundamental than sum rules derived in Ref. [161 ] . 

Our theoretical scheme is based on the BCS-Leggett 
ground state, extended [H to non-zero temperature T. 
There are then two contributions to the square of the 
pairing gap A 2 (T) = A 2 C (T) + A 2 g (T), corresponding to 
condensed (sc) and to non-condensed (pg) pairs, which 
are associated with a pseudogap. The fermions have dis- 

, — /i (In what 



persion E p = + A 2 (T), where £ p = 

follows, we omit the subscript p for convenience). 

A consistent set of diagrams which satisfy Eq.(fT|), 
include @, E3, O, both Maki Thompson (MT) and 
two Aslamazov-Larkin (AL) diagrams. In general, the 
current-current correlation function is xjj = P + ^ +W> 
where Cj is known 17, 18| and associated with excita- 
tions of the collective modes. All transport expressions 
in this paper reduce to those of strict BCS theory when 
the attraction is weak and A pg = 0. 

Consistency in linear response theory is based on the 
use of Ward identities which connect transport to the 
fermionic self energy, 



S(p,w) 



-17 - 



A 2 

pg 



cj + £p + ij cj + £p + i0 + 



(2) 



Here we use a broadened BCS self energy which was 
adopted experimentally in Refs. [19( and theoretically |6j 
to address radio frequency (RF) based studies in the cold 
gases. The condensed pairs have the usual BCS self en- 
ergy contribution, E sc , while the self energy of the non- 
condensed pairs T, pg contains an additional damping term 
parameterized by 7. 

To keep the equations simple and transparent we pro- 
ceed in two stages. We begin in the weak dissipation 
limit, by which we mean the lifetime of the fermionic ex- 
citations is very long, so that to leading order we may set 
7 rs + in Eq.©. In this weak dissipation limit, using 
the MT and AL diagrams [l7 , 18 1 one can write more 
simply 
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where K = 1, E± = E p±cl/2 , f± = f(E±) and SA 2 



A 2 - A 2 

^sc pg 



With these diagrams, one can prove consistency with 
the sum rule for the transverse susceptibility (Eq.Q). 



[For the more general case the proof depends on the fact 
that Eq.{T]) is closely tied to the absence (above T c ) and 
presence (below T c ) of a Meissner effect]. More explicitly, 
in the weak dissipation limit, we note that the total num- 
ber of particles can be written as n = J2 P (l — -^(1 — 2/)). 
The supcrfluid density at general temperatures is given 
by n s = mReP-(0, 0) + n = f £ Ep ^ (iz|£ + $L) . 
The transverse susceptibility at q — > contains no col- 
lective mode contributions; thus using P in Eq.Q, the 

left hand side of Eq.© is £ p ^ ^ 2A ™ 1 ~ 2/ 
^ 2 -AL Of' 



~W- dE 



E 2 E 

= — , thereby proving the sum 



- 4 
rule. 

We introduce the abbreviated notation 8i(u>) = 5(uj — 
E+-E-), 5 2 (u>) = 5{lo — E + + E-), and 6 is the polar 
angle. Then based on xr(w,q) the shear viscosity is 



itoj ^— v p 2 sin 2 6 



n = —m hm hm — - 

u->Oq->o 2q z ^— ' m 
p 



1 - /+ - f- 



< E+E - 7M- ~ (S^) - - (f + - /_) 



2E+E- 
E+E- 6 A 2 

2E+E_ 



(4) 



If we now take the low w, q limits in Eq.(|4|), 77 assumes a 
form similar to a stress tensor correlation function. We 
incorporate lifetime effects [l4| (which preserve the ana- 
lytic sum rule consistency) by writing 5(u ± q • V p -E) = 

After this regularization 



V 



dp 



p 6 E 2 ~A 2 pg e ( df^i 

E 2 V 



157r 2 m 2 E 2 



dE) 



(5) 



from which one can identify the effective carrier number 
(n e ff(T) oc 777 discussed in the introduction, and verify 
that it is dramatically suppressed at low T. 

Physically, the two types of terms which appear in 
Eq.Q are well known in standard BCS theory. The first 
(Si) refers to processes which require a minimal frequency 
of the order of 2A(T); they arise from the contribution of 
fermions which are effectively liberated by the breaking of 
pairs. The second of these terms, involving S 2 , arises from 
scattering of fermionic quasi-particles and is the only sur- 
viving contribution to the viscosities, which are defined 
in the u — > limit. Note that both contributions in- 
volve the difference of the condensed and non-condensed 
components (SA 2 = A 2 C — A 2 g ) with opposite overall 
signs. Importantly, the low cj quasi-particle scattering 
processes are reduced by the presence of non-condensed 
pairs - simply because more pairs lead to fewer fermions. 
By contrast in the high uj 2 A limit the number of 
contributing fermions will be increased by the presence 
of non-condensed pairs when they are broken, although 
this effect is only present in a uj ^ response. 



Figure 1: (Color online) (a) Calculated (homogeneous) viscos- 
ity r\ = nha for a unitary gas, illustrating robustness under 
changes in the value of the lifetime I/7. The color coded 
lifetimes (in units -E^ 1 ) indicated in the inset correspond to 
their counterparts in the main figure for I/7 deduced from 
fits to RF experiments (black squares) and both linear (red) 
and quadratic (blue) dependences on T/T c . (b) Contributions 
to a (thick black curve) from condensed and non-condensed 
pairs and from the simple bubble diagram. Inset plots the ef- 
fective carrier number as oc 777, showing the decrease relative 
to high T. 



Equation (J5J) is a generally familiar BCS expression 
[1, 0] except for the effects associated with non-zero A pg 

which appears as a prefactor 1 — ^f- This deviation from 
unity can be traced to the AL diagrams. Note that all 
terms are reduced by a coherence-effect-prefactor associ- 
ated with £ 2 /E 2 . The fact that the non-condensed pairs 
suppress 77 derives from the same physics discussed sur- 
rounding Eq. and associated with the scattering term 
61 . The negative sign for this bosonic term comes physi- 
cally from the fact that when pairs are present there are 
fewer fermions to contribute to the viscosity. In a more 
formal sense, this contribution is required for number 
conservation and the sum rule involving \t- 

Importantly, the expression in Eq. (|5j) can be general- 
ized to the stronger dissipation limit, by introducing gen- 
eralized Green's functions defined in Ref. ([13|), which 
represent the various pg and sc contributions, 77 — 
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pg pg 

where ± is P ± Q/2, P 
w n ,f2„ are (fermionic, bosonic) Matsubara frequencies. 
We emphasize that the lifetime, 7" 1 is a temperature 
dependent many body effect here found to be associated 
with pair-fermion inter-conversion, and partially quanti- 
fied by the cold atom Radio Frequency (RF) "photoe- 
mission" experiments 0, [ij} for 40 K. In 6 Li, where the 
viscosity experiments are performed [H, I21L 123 ] , we have 
previously estimated jl 51 ] this parameter by fitting non- 
momentum-resolved RF experiments. We note that the 
introduction of this fermionic lifetime into a broadened 
BCS form for the non-condensed pairs (or pseudogap con- 
tribution) of Eq. <j2j) is widely used in the high T c litera- 
ture p] as well. 

In Figure. QJa) we plot the shear viscosity for a unitary 



Figure 2: (Color online) (a) Comparison of shear viscosity 
77 = anh and (preliminary) experiments [13] (red triangles) at 
unitarity for a trapped gas. In theory plots (black dots) we 
use the calculated thermodynamics for the trap energy E and 
entropy density s. The inset in (a) plots data from [J], (b) 
Comparison of 77/s. The blue dashed line labels the quantum 
lower limit of 77/s given by Ref. |2j. 



gas as computed via the strong dissipation approach. We 
have verified that the results are similar if we use the sim- 
pler form of Eq.© directly. The plot is for a defined as 
rj = anh versus temperature for a homogeneous system 
at unitarity and for a range of different lifetime param- 
eterizations. The inset to Figure. Ha) presents a plot 
of the RF-deduced lifetime as black squares along with 
a few alternative functional forms. Each of these corre- 
sponds (via color coding) to the plots for a in the main 
body of the figure. In all cases a drops to zero at low 
temperatures, although one can see that this is slightly 
countered by the fact that the fermions are longer lived 
at low T. This figure should make it clear that the be- 
havior shown here is quite independent of any detailed 
models for the inverse lifetime 7 -1 . This largely follows 
from the intuition already embedded in Eq.© showing 
the dominant effect is due to exponential decrease in the 
number of condensate excitations associated with s-wave 
superfluidity. 

In the inset of Figure. |ljb) we plot the effective carrier 
number defined as oc 777 as a function of T. The curve, 
normalized to the high temperature value where A = 0, 
shows a clear suppression of the carrier number associ- 
ated with the non-condensed pairs. The main body of the 
figure presents a breakdown of the various contributions 
to the viscosity coming from this bubble term and from 
the contributions of condensed (proportional to Ag C ) and 
non-condensed (proportional to A^ 9 ) pairs. 

We turn now to calculations in a trap based on the 
local density approximation. The inset to Figure. [2^a) 
presents a plot of experimental data from Ref. [l[ The 
main body of Figure. [2{a) presents a comparison of the 
viscosity coefficient a between theory (based on the RF- 
deduced lifetime), as black dots, and experiment [22| (red 
triangles) as a function of E. Figure. [2] (b) shows the 
comparison of 77/s where s is the entropy density. We 
find that 77/s appears to be relatively T independent at 
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the lower temperatures. Here we have used the calcu- 
lated trap thermodynamics [ij to rescale the various axes, 
and our calculations are based on the same trap averag- 
ing procedure as in Ref. (22|. One can anticipate that, 
particularly at the lower T, the trap-integrated viscosity 
will be artificially higher than for the homogeneous case, 
since r\ will be dominated by unpaired fermions at the 
trap edge. Overall, it can be seen that our calculations 
agree favorably with the preliminary experimental data 
shown in the inset and figure. Interestingly, the observed 
behavior appears more consistent with previous helium-3 
experiments [l(| than those in helium-4 fllj j . 

There appear to be no other BCS-BEC based calcula- 
tions of r\ in the literature which address the entire range 
of T below T c and also the consistency check of Eq. (JTJ) . 
Bruun and Smith 1131! have studied the above T r shear 
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viscosity and recognized 13J that the pseudogap reduces 
77. However, the diagram set which was used was "not 
conserving" [13j ■ Rupak and Schafer [§] argued that r] 
is dominated by the Goldstonc bosons or phonons and 
predicted an upturn at the lowest T in both 77 and rj/s. 
This upturn has not yet been seen experimentally, and 
it seems also be missing in earlier Boltzmann-based ap- 
proaches [23j. Importantly, in BCS-based fermionic su- 
perfluids the Goldstone bosons do not couple to a trans- 
verse response such as rj. 

In summary, in this paper we have effectively as- 
cribed the low viscosity of the strongly interacting Fermi 
gases to effects associated with the depression in the 
normal fluid carrier density (n e //(T)) arising from non- 
condensed pairs and the related pseudogap. Our viscosity 
decreases monotonically as T — ¥ 0, as seen in superfluid 
helium-3 and in very recent Fermi gas data [l|, sug- 
gesting the absence of a ground state normal fluid. For 
these fermionic superfluids, lifetime effects are associated 
with fermion-boson interconversion. Moreover because of 
this mixed statistics [l2j |. Boltzmann approaches are not 
amenable to including these lifetime contributions. We 
emphasize that all approaches to transport must neces- 
sarily demonstrate, as we do here, consistency with num- 
ber conservation. 

With these on-going debates and the possible impact 
for other physics sub-disciplines, measurements and the- 
ories of the shear viscosity take on particular impor- 
tance. As shown here they constrain microscopic theories 
through the detailed information they provide about the 
nature of the superfluid excitation spectrum. Interest- 
ingly, calculations of the cj — > conductivity lead, by 
analogy, to "bad metal" behavior [24| in the cuprates, 
which may be a counterpart to perfect fluidity in the 
cold gases. 
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